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Abstract 

The celebrated Lambert-Beer law of light absorption in photochromic 
media is only valid at low intensities of incident light and low concentration 
of chromophore. Here we address the generic problem of photo-absorption 
dynamics, experimentally studying the case of azobenzene isomerization. 
We show that the non-linear regime is very common and easy to achieve 
in many practical situations, especially in thick samples where the light 
depletes the chromophore in the first layers and can propagate through 
the medium with a sub-exponential law. Importantly, the crossover into 
the non-linear absorption regime only weakly depends on the dye concen- 
tration and solution viscosity. We experimentally quantify the character- 
istics of this peculiar optical response and determine the key transition 
rate constants. 
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Introduction 

Lambert-Beer law is well known and widely used in spectroscopy: it states that 
the light propagating in a thick absorbing sample is attenuated at a constant 
rate, that is, every layer absorbs the same proportion of light [T]. This can be 
expressed in a simple form as the remaining light intensity at a depth x into the 
sample: I(x) = Iq exp(—x/D) where Jo is the incident intensity and D is a char- 
acteristic length which one calls the "penetration depth" of a given material. If 
an absorbing dye is dispersed in a solution (or in an isotropic solid matrix) this 
penetration depth is inversely proportional to the dye concentration. In this way 
it is possible to determine a dye concentration c by experimentally measuring the 
absorbance, defined as the logarithm of intensity ratio A = In (io/l) = x i c /3), 
where x is the thickness of the sample (the light path length), c the concen- 
tration of the chromophore, and 8 the universal length scale characteristic of a 
specific molecule/solvent. Thanks to Lambert-Beer law, UV-visible absorption 
spectroscopy is a useful and practical tool in many areas of science, for instance, 
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in biology to determine the degree of purity of a protein or the concentration of 
DNA [2]. One should note that in chemistry and biology one often uses base- 
10 logarithm in defining the Absorbance, Aio = \oglo/l = (ec)x. If c is in 
molar units, the constant of proportionality e is called the "molar absorption 
coefficient" and it is inversely proportional to the characteristic length S defined 
above. 

However, this empirical law has limitations, and deviations are observed 
due to aggregation phenomena or electrostatic interactions between particles. 
More importantly, in photosensitive media much stronger deviations due to 
self-induced transparency, or photobleaching [31 2] can occur. This effect has 
been reported in many biological systems like rhodopsin [SJ [5] , green fluorescent 
protein [7] or light harvesting complexes [8] . 

In most photosensitive molecules irradiation with light at a certain wave- 
length induces a conformational change (isomerization) from an equilibrium 
trans state where the benzene rings are far apart to a bent cis state where 
they are closer. The particular characteristic that makes azobenzenc and its 
derivatives interesting is that the two isomers have different absorption spectra, 
primary absorption peaks in two different regions: the trans isomer absorbs 
around 320 nm, while the cis isomer at around 440 nm [5] . The position of the 
two peaks is sensitively changed by a variety of chemical groups attached to 
the basic azobenzene molecules, making azobenzene derivatives family a model 
material for study of many photo-chromic phenomena. Irradiation with light 
at the wavelength of the trans peak progressively depletes the molecules in this 
conformation. Effective trans-cis photo-bleaching of the first layers allows a fur- 
ther propagation of light into the sample and this leads to nonlinear phenomena 
which are interesting both from the theoretical (TQl [HI El [13] and from the 
experimental point of view [TJl E3 HS1 02] ■ 

The non-Lambcrtian propagation of light through a medium has important 
consequences for the analysis of photo-isomcrization kinetics: when the photo- 
bleaching becomes important, the measured absorbance does not follow a simple 
(traditionally used) exponential law anymore. The isomerization process follows 
the first-order kinetics: 



where nt is the fraction of isomers in the trans state, kcr and kxc are the 
cis-trans and trans-cis constants of photoisomcrization, respectively, Itc an d 
Ict are the intensities of light at the wavelengths which excite the two tran- 
sitions, and 7 is the rate of spontaneous thermal cis-trans isomerization. In 
the experiments described below we use an azobenzene derivative in which the 
trans-cis and cis-trans absorption peaks are widely separated and illuminat- 
ing light monochromated at the trans-cis transition wavelength. In this case 
the stimulated cis-trans isomerization is negligible (that is, Ict ^ 0) and the 
kinetic equation reduces to 



dn t 
dt 



Irckrcnt + IcrkcT^ - n t ) + 7(1 - n t ) 



(1) 
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Figure 1: (a) Intensity I/Iq in the photo-stationary state as a function of the 
sample thickness x/D for several values of a. At small a the decay is exponen- 
tial; the light penetrates deep into the sample as a increases, (b) Intensity I/Iq 
as function of time through a fixed x/D = 2.7 and different incident intensity; as 
a grows, the dye reaches the stationary state increasingly fast with a behavior 
which differs from a simple exponential. 



In this equation the intensity / = I(x), depending on the depth into the sample. 
It is convenient to define a non-dimensional parameter a — Io^rc/l, which 
represent the balance of photo- and thermal isomerization at a given incident 
intensity Jo- I n this notation, the amount of molecules in the trans conformation 
in the photostationary state is simply 

nt = (l + al/Io)- 1 . (3) 

The change in intensity between subsequent layers of the absorbing material 
is proportional to the conversion of trans molecules in the layer. This theoretical 
problem has been studied by Statman and Janossy [12] , and Corbett and Warner 
[13] . Neglecting the stimulated cis-trans isomerization (which is appropriate in 
our study, the model can be much simplified to give: 

d ± = -l£ (4) 
dx D y ' 

where D is the penetration depth, inversely proportional to concentration. Com- 
bining (0| and © the stationary-state light intensity at a depth x is given by 
the relationship [13] 

In [I(x)/I Q ] + a[l(x)/I - 1] = -x/D. (5) 

Figurc[lja) shows the intensity variation I{x) for various values of the parameter 
a: if the incident intensity is low enough (/o^tc/7 ~~ * 0), the Lambert-Beer law 
is valid and the decrease is exponential, but if the incident intensity is high the 
bleaching of the first layers becomes progressively relevant so that they become 
partially transparent to the radiation. The decay thus tends to become linear 
in the large part of the sample, I(x) ~ /o(l — x/aD). 
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In order to model the dynamics of photoisomerization, which is evidently 
inhomogeneous across the sample, equations ([2]) and (|4]) must be coupled to 
give 



Solving this equation (see [TB] for detail) leads to the integral expression for the 
intensity I(x, t) that we would detect from the sample of thickness x at time t: 



The upper limit of this integral is the measurable absorbancc A = \yi\Iq/I{x, t)} 
from a sample of thickness x. Figure QJb) shows the resulting prediction for 
the time-evolution of intensity transmitted along the path x/D. Note that at 
t = all curves converge to the Lambert-Beer I/Iq = exp[—x/D], while at long 
times a significant portion of chromophore is bleached and transmitted intensity 
increases. 

We should note that, however important the practical case of molecular iso- 
merization might be, the problem of non-linear photo-absorption dynamics is 
much more generic. Even ordinary dye molecules that do not undergo conforma- 
tional changes stimulated by photon absorption, still follow the same dynamic 
principles, only with electronic transitions in place of trans-cis isomerization. 
Therefore, the results of this paper should be looked upon in such a generic 
sense. In particular, the two key conclusions: that the crossover intensity into 
the non-linear photo-absorption regime is independent on dye concentration 
and the rate of the transition is independent on solvent viscosity, are probably 
completely general. 

Materials and Methods 

The photochromic derivative 4'-hexyloxy-4-((acryloyloxy)hexyloxy)-azobenzene 
(abbreviated as ACgAzoCe) was synthesized in our lab by Dr. A.R. Tajbakhsh. 
The molecular structure is shown in figure [2] and the synthesis described in 
|19j . The absorption spectrum of the molecule has a peak centered at 365 nm, 
characteristic of the trans state, and a peak at 440 nm characteristic of the cis 
state. We deliberately choose this molecule with widely separated absorption 
peaks to be able to unambiguously monitor the kinetics of trans-cis transition. 



Figure 2: The monomer ACeAzoCg has the azobenzene group symmetrically 
positioned between two similar aliphatic chains, one of which is terminated by 
an acrylate group. 

The illumination was provided by a Nichia chip- type UV-LED, whose output 
power was accurately regulated by a power supply. The spectrum of emission of 
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Figure 3: I/I for 3 different values of x/D (x - 0.2, - 0.7, • - 1.1) 
and 2 different values of a, corresponding to: (a) Iq = 4/iW cm~ 2 , and (b) 
Iq = 20/iW cm~ 2 . The periodic instability was reproducible in all low- viscosity 
experiments. 



the light is centered at 365 nm (near the absorption maximum of our azobcnzcne 
derivative) and is about 10 nm wide. The LED light was attenuated by passing 
through a black tube of controlled length, placed in front of a quartz cuvette 
with 1 cm optical path. We determined the power density of light incident on the 
sample by calibration against the LED input current J [mA], the illumination 
cone angle of the LED (132°) and the illuminated area at the fixed distance 
to the sample (20 cm), giving the intensity Iq = 0.4 J /iW/cm . Several values 
of intensity were used in reported measurements, ranging between Iq = 4 and 
60 /iW/cm 2 . 

For light absorbtion measurements of a Thcrmo-Oriel MS260i spectrometer 
(focal length 260mm) was used. A shutter was placed in front of the lamp so that 
the sample can be rapidly screened from the light source. An liquid lightguide 
conducted the light from the cuvette to the spectrometer, through a 50 /im slit 
and the Andor linear-array CCD camera connected to a computer. The camera 
allowed a rate as fast as 0.021 s per frame. Before every absorption experiment, 
a background and a reference spectra were collected: the background is the 
spectrum without the illumination, and the reference is the spectrum collected 
with the LED illuminating the cuvette filled with plain solvent without the 
chromophore. The absorbancc is then calculated from the counts of the detector 
as: 

/ Reference — Background \ 



A = In 



I Signal — Background / 



All experiments were carried out at room temperature in the dark. 

For our detailed dynamic experiments, a very important issue was the vis- 
cosity of the solution. In the earlier work |19j we have shown that solvent 
viscosity has no significant effect on the trans-cis transition rate. However, at 
high illumination intensity we have encountered an unexpected problem. Figure 
[3] shows that the transmission of light through a low- viscosity dye solution (in 
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pure toluene) displays a characteristic oscillatory feature. Detailed analysis of 
this phenomenon is beyond the scope of this paper. Whether the oscillations 
are linked to the local convection due to the heating of the sample spot [5D] or 
to the diffusion of the less dense cis molecules - or whether they are intrinsic to 
the non-linear photochemical process [5T] - is not clear at this stage and would 
need further investigation. 

In order to avoid this difficulty, the dye solutions were prepared in a mix- 
ture of toluene and polystyrene of high molecular weight (Mw=280000), both 
purchased from Sigma-Aldrich and used without further purification. Adding 
polystyrene increases the viscosity of the solution by over 2 orders of magnitude, 
and in this way prevents fluid motion in the cuvette on the time scales of our 
measurements. Polystyrene-toluene solutions were prepared at a fixed weight 
ratio of 37.5wt%. 

In all cases it was important to verify that the dye concentration remained 
in a range where, at time t = 0, the linear proportionality between absorbancc 
and concentration (Lambert-Beer law) held. This is important because the con- 
centration of molecules in the trans state at every instant was determined from 
the absorbance at 365 nm. Absorbance was measured at several concentrations. 
The deviation from linearity started at A ss 3, which corresponds to the dye 
weight fraction of c = 0.03 (3 wt%) in our 1 cm cuvette. After this point, aggre- 
gation effects start playing a role and the basic Lambert-Beer law is no longer 
valid, undermining the theoretical relationship given by the equation ([7]). We 
always kept the concentrations below this value, so that the linearity at t = 
was maintained, with A = x/D where the penetration depth is inversely pro- 
portional to chromophore concentration, expressed as weight fraction, D ks 5/c 
with (5 = 91 fjja. 

We prepared three different dye solutions with (non-dimensional) weight 
fractions c = 2.5 ■ 10 -3 , 0.01 and 0.025, resulting in the values of penetration 
depth D = 36 mm, 9.1 mm and 3.6 mm, respectively (the sample thickness was 
kept constant at x = 10 mm). 

Every point of the spectrum was collected as an average of 100 measure- 
ments. All isomerization reactions were followed for several hours, until a pho- 
tostationary state was reached. The measurements were repeated at different 
illumination intensities (regulated with the power supply) and at different dye 
concentrations. 

Results and discussion 

A basic characteristic of the photoisomerization problem is the rate of spon- 
taneous thermal cis-trans isomerization 7. For a given azobenzene derivative, 
at fixed (room) temperature and sufficiently low dye concentrations to avoid 
self-interaction, this rate is approximately the same for all our solutions. We 
measured this rate after monitoring the relaxation of the spectrum after the illu- 
minating LED is switched off (see [19] for detail) and obtained 7 w 1.25'10 _4 s 
(or the corresponding relaxation time of ~ 8000s). 
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Figure 4: The effect of photo-bleaching for samples with high dye concentration 
(x/D = 2.7). Three values of irradiation intensity are labelled on the plot. Solid 
lines are fits to the data with only one free parameter a, giving a « 24.7 for the 
highest intensity, and a « 7.5 for the middle intensity. 

In order to test the predictions of the theory, dynamic absorption measure- 
ments were performed for different dye concentrations and different light inten- 
sities. Taking the equation this is equivalent to changing x/D (where D is 
inversely proportional to the dye concentration) and a, which is proportional to 
the incident intensity Jo- With this experimental setup it was possible to follow 
all the isomerization kinetics and thus the time dependence of I /Iq. 

At time t — the absorption follows the Lambert-Beer law, which guaran- 
tees the proportionality between the absorbance and the dye concentration. In 
figure 0] the representative experimental results are shown for the solution with 
the highest chromophore concentration (D = 3.6 mm, leading to x/D = 2.7) and 
three values of incident intensity. One finds that all curves converge to the same 
initial value corresponding to the I/Iq = exp[—x/D], which for this concentra- 
tion means quite a low transmission {I / Iq ~ 0.07). The classical Lambert-Beer 
would correspond to this value remaining time- independent. We fit the data 
with the theoretical model given by equation {7} where we input the values of 
7 and x/D, leaving only a = Iokxc/l free. Two data sets at higher intensity 
show the transmitted I(x, t) reach saturated value. In these case the fit is con- 
fident because we have to match both the slope and the amplitude of the curve. 
We obtain a « 24.7 for I = 60 /iW/cm 2 and a w 7.5 for I = 20 /iW/cm 2 , 
the ratio of which is close to expected 3. For the low- intensity irradiation, al- 
though the deviation from the Lambert-Beer law is evident, our fitting is not 
reliable because we did not obtain the intensity saturation in a reasonable time 
of experiment. 

At lower concentration of chromophore, corresponding to D ss 9.2 mm and 
x/D = 1.1, figure O shows the similar features of the non- linear effect. Only 
at higher irradiation intensities we could achieve the saturation and the steady- 
state value I(x) corresponding to the solution of equation ([5]). The change 
of curvature, notable in figures QJ (b) and [4] is not so clear here even at the 
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Figure 5: The same experiment as in figure [4] but with an intermediate dye 
concentration (x/D = 1.1), and the same values of irradiation intensity. Here 
the solid lines are not fits, but theoretical plots of equation ([7]) for a = 24.7, 7.5 
and 1.64 for the decreasing intensity, respectively. 



highest Iq. However, in the comparative analysis of data we now take a different 
approach. Assuming all the parameters for the curves are now known (7 and 
x/D from independent measurements, and a from the fitting in figure 0]), we 
simply plot the theoretical equation ([7]) on top of the experimental data. For 
the low-intensity curve we take a ~ 1.64, which is a 1/5 fraction of the value for 
the highest Iq assuming the linearity. It is clear that the theory is in excellent 
agreement with the data. 

Finally, we examine the dynamics of non-linear absorption at low dye con- 
centration (D ss 91 mm, x/D = 0.27) in figure [6j In this case the complicated 
integral equation ([7]) simplifies dramatically [18], since at small i/D < 1 the 
difference between A = ln[I /I] and x/D (which is the range of integration in 
([7]), is also small. The integration can then be carried out analytically, giving 



In 



I(x,t)\ x 



1 - - 1-e 



-7(l+a)t 



(8) 



which gives in the stationary state the correct solution of equation (J5j) approxi- 
mated at small x/D: 



hi 



I(x)\ x/D 



Io J 1 



More importantly, we also see that that the rate of the process described by 
the approximation (JSJ) is given by the simple exponential, t -1 = 7(1 + a) = 
7 + kxclo- This is in fact the rate originally seen in the kinetic equation @. 
Therefore, we now fit the family of experimental curves in figure [5] by the sim- 
ple exponential growth and then analyze the relaxation rates obtained by this 
fit. The inset in figure [6] plots these rates for all the io values we have stud- 
ied. A clear linear relation between the relaxation rate and Iq allows us to 
independently determine the molecular constant: kxc ~ 10 _4 cm 2 s _1 /uW . 
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Figure 6: At low dye concentration (x/D = 0.27) the sample is relatively trans- 
parent. The data are for the same three values of irradiation intensity as in 
the earlier plots. The inset shows the plot of exponential relaxation rate t -1 
against Iq, with the linear fit. 

Conclusions 

The main conclusion of this paper is that one has to be cautious with the classical 
concept of light absorption, represented by the Lambert-Beer law. The ques- 
tion is not about more complicated multi-particle effects at high concentration 
of photo-sensitive molecules in the system: obviously for very high absorbancc 
values one does not expect linearity. Importantly, even at very low concentra- 
tions (corresponding in our case to the low x/D ratio) the illumination intensity 
above a certain crossover level would always produce a non-linear dynamical 
effect equivalent to the effective photo-bleaching, increasing the effective trans- 
mittance of the sample. The crossover is expressed by the non-dimensional 
parameter a = Iokrc/l > 1 and is, therefore, an intrinsic material parameter 
of every chromophore molecule, but not dependent on the dye concentration. 
Note that the thermal cis-trans isomcrization rate 7 is strongly temperature 
dependent, influencing the crossover intensity. 

The experiments, carried out to quantify this dynamical effect, have con- 
firmed the model predictions and demonstrate several linked features, especially 
pronounced in the high-concentration systems. The final saturation steady-state 
value of transmitted intensity (or the practically measured absorbance) is linked 
with the rate at which the steady state is reached. At high intensity of illumi- 
nation one finds a characteristic sigmoidal shape, with the changing curvature 
of I(t) variation. Our experiments were deliberately carried out in a highly 
viscous solvent to eliminate the additional complexities caused by the possible 
local convection flows of different isomers, which arc illustrated in figure [3] Cer- 
tainly, a much more in-depth study will be required to take such effects into 
account. 

It is important to emphasize that these results and conclusions are very 
general, not at all confined to the class of azobenzene or other photo-isomcrizing 
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molecules. All dyes are essentially systems with two long-lived energy levels and 
thus the non-linear kinetics of their absorption and saturation will obey the same 
rules as discussed in this paper. 
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